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ABSTRACT 

In this study, turbulent flow in a curved channel of 
square cross section was investigated experimentally to 
determine secondary flow effects. Probe surveys were con- 
ducted to establish vortex and total pressure loss distri- 
bution in the exit plane of a 90^ bend. Overall losses 
were determined by measuring the momentum of the flow with 
a force plate for different Reynolds numbers at turning 
angles of 0°, 45°, 90°, and 135°. 
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INTRODUCTION 



I • 

The art of designing tur bomachines has improved radi- 
cally during the past decades, A better understanding of 
flow conditions in rotors and stators has been responsible 
for much of this improvement. Continued advancement in 
turbomachinery cannot be achieved without an increased 
knowledge of flow patterns and conditions which exist in 
various blade configurations. 

An important flow phenomenon known to exist but not 
too well understood is that of secondary flow produced 
when flow is directed between the curved surfaces of adja- 
cent blades in a cascade. Losses caused by secondary 
flow reduce the efficiency of turbomachines, A better 
understanding of these losses may result in establishing 
design techniques to minimize them. 

The purpose of this study was to investigate secondary 
flow effects produced when a turbulent flow is directed 
through a curved channel of square cross section. As 
shown in Figure 1, the vertical walls of the channel repre- 
sent the surfaces of adjacent blades in a cascade. Pro- 
cedures required to investigate secondary flow in a 
cascade test rig are complicated, and the distinguishing 
of secondary flow effects from tip effects, trailing vor- 
tices, and wakes can be very difficult. The simple model 
of Figure 1 eliminates many of these complications. 
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Probe surveys were conducted to investigate the flow 
field at the exit of a 90^ bend. Overall secondary flow 
losses were determined for bends with different turning 
angles by using a force plate to measure the useful exit 
momentum. To accomplish this, a hydraulic force plate 
apparatus was designed. 
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Representation of Adjacent Blade Surfaces by 
a Square Channel 
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APPROACH TO PROBLEM 



II . 

A. SECONDARY FLOW EFFECTS 

The results of applying the equation of motion for 
incompressible flow 

V(|-) = V X (V X V) 

to flow passing through a rectangular curved channel is 
taken from Vavra [Ref. 1] and shown in Figure 2. The vor- 
tices represent the well known secondary flow. Losses 
associated with secondary flow can be separated into two 
categories . 

Flow in a straight channel has total pressure losses. 
For the flow velocities and channel size being considered, 
these losses tend to be confined to the boundary layer 
regions along the walls and in the corners. In curved 
channels, secondary flow vortices distribute boundary 
layer flow throughout the channel. This results in a 
higher total pressure loss for a curved channel than for 
a straight channel. 

The second type of loss is due to angular deflection 
of the flow in different directions. The direction of the 
overall path of the flow is obtained by averaging individ- 
ual particle velocities over the channel cross section. 
Velocity components perpendicular to the overall flow path 
are produced by secondary flow. In the blading of turbo- 
machines, these velocity components do not contribute to 
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Figure 2 

Secondary Flow Vortices in a R.ectangular Bend 
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the energy conversion in the machine and can be classified 
as non-usef ul • 

Analytical treatments of secondary flows in turbo- 
machines were done in Ref. 3. It was concluded that no 
method was available to adequately predict losses caused 
by secondary flow. It is desirable^ therefore, to investi- 
gate and measure these losses in well controlled experi- 
ments . 

B. FORCE PLATE DETERMINATION OF LOSSES 

An isentropic process from uniform conditions at sta- 
tion 1 to station 2 of Figure 2 would yield an exit veloc- 
ity of . A process with total pressure losses between 

"^th 

1 and 2 would yield an average velocity > which is 

smaller than . The difference between and in- 

th ^ ^th 

creases with increasing loss. The velocity distribution 

at station 2 can be determined by probe surveys. A com- 
plete survey would be required for each turning angle and 
flow rate investigated. The impact force of the flow at 
station 2 is directly related to its average velocity. A 
large plate positioned perpendicular to the flow path as 
shown in Figure 3 serves to measure the impact force of 
the flow. Such a plate measures only the velocity compo- 
nents parallel to the flow path, thus eliminating the non- 
useful velocities. Conservation of momentum is applied to 
the control volume shown in Figure 3. The flov; enters the 
control volume through surface 2 and is assumed to exit 
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Figure 3 

Conservation of Momentum Applied to the 
Force Plate 
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perpendicularly through control surfaces 4 and 5* The 
static pressure along surface 6 is assumed to be atmo- 
spheric (P ) . The force on the control volume in the 
x-direction must be equal to the change of momentum in 
the x-direction. Thus 



- \ pV^dA = \ P„dA - \ P_dA + \ P dA 

Ja ^ Ja ^ ‘^A ^ '^A ^ 

^2 2 3 6 



The force on the plate, F, is given by 



=1 



PjdA - 



( 1 ) 



( 2 ) 



or 



I 



P,dA = F + P,A„ 
3 A3 



Substituting this relation into (1) 



[ 1 ‘'2“ - + V3’ + Vc 



- 1 = S. 



"^2 ^2 



F = \ pV^dA + V P dA - P.A 



F = \ pV^dA + \ (P^- P.) dA 

J X X ^ 2 A 

X 2 



(3) 
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Probe and pressure tap measurements indicated that was 

equal to P^ within 0,06 inches of water in all experiments 



using the force plate. Based on this experience, P^ was 



A loss coefficient is defined in Appendix A using force 
plate measurements . 

C, PROBE SURVEY OF THE FLOW FIELD 

Determination of secondary losses is important in iden- 
tifying their magnitude in channels with different turning 
angles and Reynolds numbers, A better understanding of 
the secondary flow phenomena may be achieved if these mea- 
surements are combined with probe surveys of the flow field 
at the channel exit. The secondary flow vortices shown in 
Figure 2 are the result of theoretical analysis. The actual 
location of vortices can be determined only by experiment. 



assumed to be equal to P • Equation (3) then becomes 

Ix 




15 





« 



•J 




Ill* DESIGN OF FORCE PLATE APPARATUS 



A. PRELIMINARY CONSIDERATIONS 

An apparatus was needed to measure the force exerted 
on the force plate by the flow. A hydraulic apparatus was 
devised for this purpose which is shown in Figure 4. The 
force plate is directly connected to the piston and the 
oil in the cylinder is assumed to be incompressible. Any 
change in the force F will result in an instantaneous 
change in P^ . The leakage at represents the piston- 

cylinder leakage and is a function of P^^ • If P^^ increases, 

leakage through A^ increases while flow through A^^ decreases. 
Hence the volume of oil in the cylinder decreases, causing 
the piston to move in the y-direction. As the piston moves, 
the valve at A^^ opens, admitting more oil to the cylinder 
to compensate for the higher leakage. The interactions are 
the same but reversed if F decreases. Three considerations 
which lead to the, final design are discussed in the sections 
which follow. 

1 • Stability 

Before establishing a mechanical design, the system 
of Figure 4 was investigated for stability. At any instant, 
the flow rate through A^ must equal the flow rate through 
plus the rate of increase in system volume. Hence 

g + yH(t)=Q(t) 
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Figure 4 

Operating Principle of Hydraulic Force Plate Apparatus 
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where 




^0 = F 



when t = 0 



Equation (5) has the solution 

t t 

- V H(t)dt *\ H(t)dt 

Jq p c Jq 

y(t)= e \ Q(t) e 



dt + C e 



H(t)dt 

"^0 



( 6 ) 



where 

" = ^0 

System response was evaluated for the following conditions 
and forced pressure variations. Solutions were obtained 
by numerical trapezoidal integration. 

= 0.203 Ib/in^ for t < 0 

= Pq + A Pq sin (2 tt ^ i ® 

y^ = 0.001 in 

Pq = 5.0 Ib/in^ 

X = 0.5 

w = 1 . 0 in 

A = Att in^ 

P 

3 

p = 0.921 gm/cm (Turbine Oil Esso 15) 

P = 200 Ib/in^ 

P 

= 0.0314 in^ 
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Case I: T = 0.05 sec 

Case II: T = 0.10 sec 

Case III: T = 0.20 sec 

The results are shown in Figures 5-7. The disturbances 
at 0.25 seconds in Figure 5 and 0.07 seconds in Figure 6 
were attributed to the numerical method used to solve 
Equation 6. From these results, the apparatus appeared 
to be stable. 

2 . Friction 

After it was verified that the system was stable, 
a mechanical design was undertaken. The initial proposal 
is shown in Figure 8. Motion of the piston due to a change 
in F changes the amount of oil admitted from the pump. For 

each force applied, the piston moves to an equilibrium posi 

tion such that the flow entering the cylinder equals the 
leakage flow. The stopper arrangement prevents the plate 
and piston from rotating. 

Friction between the piston and cylinder decreases 
the measuring accuracy. Influences of static friction be- 
tween piston and cylinder surfaces were eliminated by rotat 
ing the cylinder. The design of Figure 8 was modified to 
allow for cylinder rotation by a low-speed electric motor 
and is shown in Figure 9. The addition of oil grooves and 
outer housing were necessary to allow for cylinder rotation 
Although the design of Figure 9 eliminated static friction, 
it increased the problem of leakage. 
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Figure 5. Piston motion "y" as a function of time "t" for Case 
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Piston motion for Case II 
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Piston motion for Case III 




Figure 8 

Original Design Configuration of the Force Plate Apparatus 
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Force Plate Apparatus Modified for Cylinder Rotation 
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3 • Leakage 



Leakage between the piston and cylinder was neces- 
sary to provide proper lubrication in minimizing friction. 
Additional leakage between the cylinder and housing intro- 
duced by the design of Figure 9 was eliminated by 0-ring 
seals as shown in Figure 10, Upon completion of the force 
plate apparatus, problems were encountered with proper 
0-ring sealing and are discussed in Appendix B, 

B. FINAL DESIGN 

These considerations, along with other factors, led to 
the final design of the force plate apparatus. The design 
drawings and photographs of the finished parts are shown 
in Figures 11 - 17. All parts were made at the Naval Post- 
graduate School, Monterey, California. Grades of steel 
specified on the drawings were not used. The piston was 
made of stainless steel whereas other parts were made of 
steel whose grades are not known exactly, to avoid delays 
associated with ordering specific materials. 

The oil catcher photograph (Figure 15B) shows an addi- 
tion hot shown in the drawing. This device was welded in 
place after machining and connected to a suction pump for 
drainage purposes. The grease fittings in Figure 14A pro- 
vided lubrication for areas not lubricated by oil. 
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Position of 0-ring Seals 
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Figure 11 
Assembly Drawing 
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Piston Drawing 




Figure 12B 
Piston Photograph 
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Figure 13A 
Cylinder Drawing 





Figure 13 B 
Cylinder Photograph 
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Figure 14B 
Housing Photograph 
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Figure 15A 
Oil Catcher Drawing 




Figure 15B 

Oil Catcher Photograph 
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Figure 16A 
Back Plate Drawing 




Figure 16B 

Back Plate Photograph 
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Force Plate Drawing 
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IV • EXPERIMENTAL APPARATUS, PROCEDURES, AKD RESULTS 

The experimental set-up is shown in Figure 18, Air 
supplied by the All i s -Chalmer s twelve-stage axial compres- 
sor passed through a control valve and into a settling 
tank. The flow then passed through a standard ASME sharp 
edge orifice into a second tank. From this tank, the flow 
passed through an eight-inch diameter pipe into a transi- 
tion piece which changed the cross section to a five-inch 
square area. From this point, the flow passed through a 
honeycomb flow straightener and into the model. Figure 18 
shows the straight channel model configuration. All chan- 
nels of different turning angles had a centerline radius 
of 12,5 inches and were attached at the connection point. 
The flow rate was controlled by the valve and determined 
from measurements taken with two slant tube mercury mano- 
meters combined with a total temperature measurement at the 
orifice. Total pressure and temperature were measured in 
the eight-inch diameter pipe. All pressures measured be- 
tween the second settling tank and the model exit were 
made with water U-tube manometers. The mercury manometers 
had an accuracy of + .005 inches of mercury, U-tube mano- 
meters had an accuracy of + ,05 inches of water, 

A. PROBE SURVEYS 

The first set of experiments consisted of Prandtl 
probe surveys at the discharge of a 19,5 inch long straight 
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Experimental Set-Up 





I. 








channel to investigate the velocity profile and to deter- 
mine the mass flow rate* All channel lengths were measured 
from the connection point to the exit along the channel 
centerline* With the exception of regions near the wall 
and in the corners, the velocity profile was nearly uniform 
as shown in Figure 19. The mass flow obtained by numerical 
integration of the probe measurements over the channel exit 
agreed with the orifice measurement within 1.5 percent. 

The second set of experiments were performed with a 90^ 
bend approximately 20 inches long. The purpose of these ex- 
periments was to investigate the flow field existing at the 
bend discharge. Yaw angles were determined with a Cobra 
probe and are shown in Figure 20. The arrows represent the 
angular deflection of the flow with respect to a plane per- 
pendicular to the exit plane. Survey points along the lower 
wall are absent because of probe vibrations. The probe was 
reinforced to prevent this problem in later experiments. 

The results showed an overall angular deflection in the x- 
direction. This indicated that the flow was not turned by 
90^ but was following the path shown in Figure 2.1. As a re- 
sult, secondary flow vortices were imposed upon the average 
angle 0. To eliminate this condition, a straight section 
10 inches in length was attached to the discharge of the 
90° bend. 

The third set of experiments was carried out with this 
arrangement. Yaw and pitch angles were determined with a 
Cobra probe mounted in the horizontal and vertical positions. 
The results are shown in Figures 22 and 23. The arrows 
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Figure 19 

Prandtl Probe Survey of the Straight Channel 
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Figure 20 

Yaw Angles in the 90 Degree Bend 
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Figure 22 

Yaw Angles Induced by Secondary Flow in the 90 Degree Bend 
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Figure 23 

Pitch Angles Induced by Secondary Flow in the 90 Degree Bend 
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represent flow angles measured with respect to a plane per- 
pendicular to the exit plane. From these results, the 
angle 0 of Figure 21 appears to be equal to zero. 

The same arrangement was used for the final set of sur- 
veys which were taken with a Kiel probe. The purpose of 
these experiments was to determine the total pressure loss 
distribution over the exit plane. The loss is represented 
by the factor 




where = reference total pressure measured by a Kiel probe 
in the eight-inch diameter pipe, 

= total pressure measured by the survey Kiel probe 
at the model exit. 

The magnitude of X does not represent total pressure loss 
caused by secondary flow. Losses that occur between the 
reference Kiel probe and the connection point are also in- 
cluded in X along with secondary losses from the connection 
point to the channel exit. Values of X only important 

in showing the distribution of the losses over the exit plane. 
The results are shown in Figure 24, The curves represent 
contours of constant loss, 

B. FORCE PLATE MEASUREMEKTS 

The hydraulic force plate apparatus designed for this 
study was not used. Problems encountered with the apparatus 
are discussed in Appendix B, The force plate system that 
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Contours of constant total pressure loss with 
indicated . 



Figure 24 

Total Pressure Loss Distribution for the 90 Degree Bend 
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was used is shown in Figure 25% The force on the plate 
was transmitted to a variable reluctance force transducer 
by a steel shaft supported by axially sliding ball bear- 
ings, The force transducer was loaded in tension and 
calibrated to display force in pounds and tenths of pounds 
on the counter. Adjusting bolts were used to adjust the 
height and pitch of the plate. The maximum system error 
was + 0.2 pounds . 

Experiments were conducted with channels of 0^, 45^, 

90^ and 135^ of turning. All channels were approximately 
thirty inches in length and the last ten inches of each 
channel consisted of a straight section. Investigation of 
the flow at the discharge indicated that the overall flov7 
path was perpendicular to the exit plane. The force plate 
was positioned thirteen inches from the exit plane and 
parallel to it. Loss coefficients were determined as de- 
scribed in Appendix A. The loss coefficients are tabulated 
in Table I and plotted against Reynolds number in Figures 
26 - 29 with the Mach number for each point indicated, ^ 
is the average total loss coefficient for a particular 
channel. Average secondary flow loss coefficients (Cg) for 
the curved channels were taken as the difference between 
the average total loss coefficient for each bend and the 
average total loss coefficient for the straight channel. 
Hence, (45°) = 0.0166 

r (90°) = 0.0317 

o 

(135°)= 0.0722 
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Force Plate Apparatus 
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Figure 27 

Loss Coefficients for the 45 Degree Bend 
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Figure 28 

Loss Coefficients for the 90 Degree Bend 
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Loss Coefficients for the 135 Degree Bend 



CONCLUSIONS 



V* 

The force plate was extremely valuable in determining 
the actual momentum of the flow. The measurements were 
very successful in determining average loss coefficients. 

As shown on page 50, there is a definite increase in the 
average total loss coefficient with turning angle. No 
assumptions were made concerning velocity profiles, compres- 
sibility, or adiabatic flow. These and other assumptions 
were only necessary to determine the theoretical momentum. 

A definite dependence of loss coefficient on Reynolds num- 
ber and Mach number cannot be established from the results 
obtained. Force plate msasurements over a wider range of 
flow rates may yield such a dependence. 

The probe survey results of the 90° bend shown in Fig- 
ures 22 - 24 raise questions which can only be answered by 
further investigation. Angular deflections in the lower 
right corner of Figures 22 and 23 are opposite to the 
theoretical vortices in Figure 2. The location of the low 
total pressure loss region in Figure 24 along with a com- 
plete absence of loss concentrations in the right upper 
quarter of the channel seem peculiar. These results may be 
due to changes in the flow field caused by the probe. Yaw 
angles and total pressure loss distributions were determined 
with probes mounted on the upper wall while pitch angles 
were determined by a probe mounted on the left wall of 
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Figure 23. Surveys taken by probes mounted on the lower 
and right walls may yield different results. 
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APPENDIX A 



REDUCTION OF FORCE PLATE MEASUREMENTS 



Applying conservation of momentum to the force plate 
resulted in Equation (4), or: 



F = - pV^dA 

g J. X 



where 



g = 32.174 



f t -Ibm 



lb f -sec 

F = force on the plate in Ibf 

p = density at model exit in ^ 

ft 



ft 



V = axial velocity at model exit in 

X S G C 

A = cross sectional area at model exit 
2 

in ft 



F represents the integration of momentum over the exit 
plane of the model and is used to determine the following 
loss coefficient (C) • 






pV^dA 

X 



£ = 1 - 



1 ‘■''th" 



^th exit velocity that could be achieved 

from the isentropic expansion shown in Figure 29% Station 
1 is shown in Figure 18 and station 2 is at the model exit. 
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T 

T1 



T 



th 




Figure 30 

Entropy Diagram for an Isentropic Process 



th 

J 

C 



= total temperature at Station 1 in R 

Ibf 



= total pressure at Station 1 in 



ft 



= atmospheric pressure in 



Ibf 

ft^ 



= theoretical static temperature at Station 2 



in °R 



= 778.16 



ft-lbf 

Btu 



= 0.24 



Btu 

lbm-°R 



T was assumed to be equal to the total temperature mea- 
1 

sured in the circular pipe. Static pressure at station 2 
was assumed to be uniform over the exit plane and equal 
to atmospheric pressure. Conditions at station 1 were 
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assumed to be uniform over the channel cross section as 
were theoretical conditions at station 2. was deter- 

mined from the average value of four wall static pressure 
tap measurements taken at station 1* With these assump- 
tions and measurements* \ pV^, dA was determined as 

j th 

follows : 



where 



”■ ■ " RIj 



R = 53.34 



f t-lbf 
lbm-°R 



25 2 

A = A = ~ — f t ^ 
^2 144 “ 



and 



V 



T = T 
1 T, 



Since 



m = 



PlA^Vi 



R(T - V^/B) 
1 



f ''i 



and therefore 









2 4*2 

+ I (mR)^ 



2mR 



The isentropic relationship gives 



P = P 
T 1 

il 



7-1 
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«£■ 



I 



t: 



I 







where it is assumed that 



7 



1.4 



Then , 



where , 



or 



Now , 






pV^^dA 






V 



2 

2 

th 



= T - T 

B th 






7-1 

7 n 



1 - 



so that , 



th RT 



th 



RT, 



. [ y 



7-1 

7 



7-1 



Pth ^h ^2 = 




The final expression for ^ is therefore: 

7-1 



? = 1 - 



'• '[y 



7 “1 
7 



■ a) ] 



B 1 



/ I ^A^2 
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APPENDIX B 



PROBLEMS ENCOUNTERED WITH THE FORCE PLATE APPARATUS 

AND RECOMMENDATION 

The first problem that occurred with the hydraulic 
force plate apparatus was excessive oil leakage% The 
0-ring grooves were designed for a rotary 0-ring seal. 
Teflon 0-rings ware chosen to keep friction betvjeen the 
cylinder and housing as low as possible. The 0-rings were 
ordered from several companies but delivery times were 
excessive. Rubber 0-rings were obtained but failed to 
provide a good seal. To help solve this problem, a new 
cylinder was made. The machinist was instructed to omit 
the 0-ring groove and make the piston-cylinder and cylin- 
der-housing clearances as small as possible without 
restricting rotation. During machining, the inner surface 
of the housing was found to be slightly out of round. This 
limited the housing-cylinder clearance while a very close 
fit was obtained between the piston and cylinder. The 
apparatus was assembled with an 0-ring placed in the hous- 
ing 0-ring groove. An excellent seal was obtained due to 
the decreased cylinder-housing clearance. The apparatus 
was tested without the force plate attached and the suction 
pump was able to handle the oil leakage. 

The force plate was attached to the piston and the oil 
pump and electric motor were turned on. Within a few 
seconds, the piston was tightly frozen in the cylinder. 
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The apparatus was disassembled and the piston-cylinder 
surfaces were badly scored. Apparently, a metal particle 
became lodged between the two surfaces during rotation. 

The oil supplied from the pump passed through a filter 
before entering the apparatus although the type of filter 
and its condition were not known. The weight of the 
force plate (ten pounds) produces areas where the piston- 
cylinder surfaces are forced together. These locations 
are especially susceptible to scoring due to surface rough- 
ness or particles in the oil. The following recommenda- 
tions are made for further work with the hydraulic force 
plate apparatus: 

1, An investigation be made to determine the type of 
oil filter system needed, 

2, The out-of-round housing be corrected, 

3. Clearances between the second cylinder, housing, 
and piston be maintained, 

4. The cylinder 0-ring groove in Figure 13A be added 
to the second cylinder. 

5. The rubber 0-rings (Porter Seal 568-040 and 568-041, 
Buna-hardness 70) be used, 

6, The piston, cylinder, housing, and back plate be 
made from materials specified on the original drawings, 

7, The force plate be re-designed for minimum weight, 

8. The piston-cylinder surfaces be checked for a proper 
finish especially where holes have been drilled. 
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